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ABSTRACT 

Context. While the halo mass function is theoretically a very sensitive measure of cosmological models, masses of dark-matter halos 
are poorly defined, global, and unobservable quantities. 

Aims. We argue that local, observable quantities such as the X-ray temperatures of galaxy clusters can be directly compared to theo- 
retical predictions without invoking masses. We derive the X-ray temperature function directly from the statistics of Gaussian random 
fluctuations in the gravitational potential. 

Methods. We derive the abundance of potential minima constrained by the requirement that they belong to linearly collapsed struc- 
tures. We then use the spherical-collapse model to relate linear to non-linear perturbations, and the virial theorem to convert potential 
depths to temperatures. No reference is made to mass or other global quantities in the derivation. 

Results. Applying a proper high-pass filter that removes large enough modes from the gravitational potential, we derive an X-ray 
temperature function that agrees very well with the classical Press-Schechter approach on relevant temperature scales, but avoids the 
necessity of measuring masses. 

Conclusions. This first study shows that and how an X-ray temperature function of galaxy clusters can be analytically derived, avoid- 
ing the introduction of poorly defined global quantities such as halo masses. This approach will be useful for reducing scatter in 
observed cluster distributions and thus in cosmological conclusions drawn from them. 

Key words, cosmology: theory - methods: analytical - galaxies: clusters: general - dark matter - X-rays: galaxies: clusters - cosmo- 
logical parameters 



1. Introduction 

Populating the far end of the halo mass function, galaxy clus- 
ters are in principle highly sensitive indicators of the cosmo- 
logical parameters and non-linear structure growth. Combining 
Gaussian random density fields with linear structure growth and 
spherical collapse, the Press-Schechter mass function and its 
variants turn out to reproduce the halo mass function in fully 
non-linear cosmological simulations extremely well. If measur- 
able, the abundance of halos in the exponential tail of the mass 
function and its evolution on cosmic time scales allow precise 
constraints on both the density-fluctuation amplitude today and 
during the second half of the cosmic age and on the matter- 
density parameter. The exponential dependence of the abun- 
dance of massive halos on cosmological assumptions promises 
tight constraints. 

A direct comparison between the theoretically predicted 
mass function of massive halos and the observed distributions 
of galaxy clusters in various observable quantities, such as the 
flux and the temperature of their X-ray emission or the veloc- 
ity dispersion of their member galaxies requires observables to 
be translated into mass. While this conversion appears straight- 
forward under the idealised assumptions of spherical symmetry, 
thermal, and hydrostatic equilibrium, the cluster population as a 
whole shows all signs of being dynamically active. It is doubtful 
whether precise cosmological conclusions can be drawn based 
on symmetry assumptions. 

Even if clusters satisfied the idealising assumptions typically 
underlying their cosmological interpretation, their mass is not an 



observable. In fact, the mass of a dark-matter halo is a poorly 
defined, derived quantity to which hardly any precise meaning 
can be given. It is common to operationally define halo masses 
as enclosed by spheres containing an average fixed overden- 
sity. However, the many different choices of apparently appro- 
priate overdensity values in the literature demonstrate that there 
is no uniquely defendable choice. If the overdensity is chosen 
very high, the masses obtained are core masses rather than halo 
masses, and if it is chosen low, density profiles constrained near 
the core need to be extrapolated into regions where they are typ- 
ically poorly measured or not at all. 

Halo definitions in numerical simulations illustrate the same 
problem in a different way. There, halos are typically identi- 
fied by group finders connecting particles with neighbours closer 
than a certain linking length. Recipes exist for how the linking 
length should be chosen, but there is no objective criterion. The 
dependence on the linking length may be less relevant in practice 
because halo masses can again be defined as the masses of all 
particles in spheres containing a fixed overdensity. However, this 
refers back to the largely arbitrary overdensity threshold and cre- 
ates the additional problem that several different plausible defi- 
nitions of halo centres exist that often yield discrepant results. 

Three main classes of observation are used to constrain clus- 
ter masses: gravitational lensing, X-ray flux and temperature, 
and galaxy kinematics. None of them measures cluster masses. 
Gravitational lensing measures the curvature of the projected 
gravitational potential. X-ray observables are primarily deter- 
mined by the gas density and temperature, which respond to 
the depth of the gravitational potential and its gradient, as do 
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galaxy kinematics. Thus, cluster observables constrain the grav- 
itational potential rather than any kind of mass. The conversion 
of the potential into a mass is hampered by the fact that mass 
is a non-local quantity, requiring an integration over potential 
derivatives. We raise the question whether cosmological conclu- 
sions can be drawn directly from cluster observables without the 
detour through problematic definitions of cluster masses. 

As one step towards a possible answer, we derive here the X- 
ray temperature function from a locally defined quantity, namely 
the gravitational potential. To this purpose, we first derive a func- 
tion predicting the number density of potential minima having a 
certain depth. We include the non-linear evolution of the poten- 
tial by considering the collapse of a spherical and homogeneous 
overdensity, and locally relate the non-linear potential depth to a 
temperature using the virial theorem. This direct relation of the 
temperature to the gravitational potential allows us to avoid in- 
troducing a global quantity such as the mass and the ambiguities 
in its definition. The formalism proposed in this work may thus 
contribute to reducing the systematic uncertainty in comparisons 
between theory and observation by avoiding empirical relations 
between cluster masses and observables. 

Unless declared otherwise, we shall use the following cos- 
mological models and parameters: Einstein-de Sitter (EdS) 
Q m0 = 1.0, Q A o = 0.0, Q b o = 0.04, h = 0.7, cr 8 = 0.52; ACDM 
D. m0 = 0.3, Q A0 = 0.7, Q b0 = 0.04, h = 0.7, <x 8 = 0.93; OCDM 
Q m0 = 0.3, Q A0 = 0.0, Clbo = 0.04, h = 0.7, <x 8 = 0.87. The 
different values for erg reflect the normalisation of the power 
spectrum to rep roduce the local abundance of galaxy clusters 
(lEkeet al.lll996h . 

2. Gaussian random fields 

Simple models of inflation predict that the density contrast, 
6(r) = p(r)/pb - 1, where p(r) is the actual density at posi- 
tion r and pb is the mean cosmic background density, should 
be a Gaussian random field right after inflation. Since the den- 
sity contrast and the gravitational potential are linearly related, 
the latter is then also a Gaussian random field. 

In this section , we sh all follow the line of argument presented 
by iBardeen et alj (Il986l) . For better comparison to their paper, 
we adopt F{r) for the random field and r)(r) = VF(r) and £y(r) = 
didjF(r) for its first and second derivatives, respectively. 

2.1. Definition 

An n- dimensional random field F(r) assigns a set of random 
numbers to each point in n-dimensional space. A joint proba- 
bility function can be declared for m arbitrary points r, as the 
probability that the field F, considered at the points rj, has val- 
ues between F(rj) and F(rj) + dF{rj) with j = l,...,m. 

A Gaussian random field is a field whose joint probability 
functions are multivariate Gaussians. Let y, with i = 1 , . . . , p 
be a set of Gaussian random variables with means (y,-) and 
Ay,- := y,— (yi). The covariance matrix M has the elements My : = 
(AyiAyj), and the joint probability function of the Gaussian ran- 
dom variables is 



P(y l ,...,y p )dyi---dy p 



with the quadratic form 



1 



VW det(M) 



e: =l£Ay,(M->) Ay,. 



e- Q d yi ---dy p (1) 



(2) 



A homogeneous Gaussian random field with zero mean is fully 
characterised by its two-point correlation function %{r\,r2) - 
%(\r\ - '■- {F{_r\)F{r2}) or equivalently its Fourier transform, 
the power spectrum P{k). 

2.2. The minimum constraint 

An expression for the number density of field minima can be de- 
rived as follows. The joint probability function for a Gaussian 
random field in three-dimensional space with zero mean includ- 
ing first and second derivatives reads 



p(F,t],OdFd 3 T]d b ( = 



1 



: e- Q dFd 3 r]d b ( , (3) 



V(27r) 10 det(M) 

with the quadratic form Q given in Eq. © and y = 

(F, 771, 772, 773, £11, £22. &3. £12, £13, 

The matrix M contains all auto- and cross-correlations be- 
tween these quantities, which read 

~-2 



(FF) = 0% , 


(Vflj) = 






TO = -y *y . 


<&&> = 


{FVij) = , 





1 



15 
. 



(4) 



where the cri, < j < 2, are the spectral moments of the power 
spectrum P(k), 



(5) 



The Fourier transform of the top-hat window function 

with the filtering scale R is 

3 (sin u - u cos u) 



W R (k) = 



with u-kR. 



(6) 



Let ro be a minimum of the field F, hence 7/(ro) = 0, 
and the eigenvalues (£1, £3) of the tensor of second 
derivatives must be positive. Within an infinitesimal volume d 3 r 
around r — 0, we can approximate 77, m gyrj and thus replace 
d 3 ?7 = I det(^)|d 3 r in Eq. (|3). We also transform the volume el- 
ement in the space of second derivatives, d 6 £, into the space of 
eigenvalues, 



d 6 ^- = t - h) (h - h) (?i - h)\ ^hAh 



(7) 



(IBardeen et al.lll986l) . Using det(£) = l\hh, we a 11 "^ at the 
final equation for the number density of minima, 

n(F) = ^ Jd^J d^ 2 J dh llr&U (8) 


x |(?i - h) (h - h) (?i -h)\p (f, r, = ^h,li, h) ■ 

2.3. Number density of potential minima 

We now apply this formalism to the Gaussian random field of 
gravitational-potential fluctuations <t>. We continue using rj for 
the first derivative of the field and f for its tensor of second 
derivatives, but introduce new variables. Instead of the eigen- 
values Xi with i - L 2, 3, we switch to the linear combinations 



Ad) + &+£» 
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These choices simplify the correlation matrix M, and we can 
later easily identify the Laplacian of the field. In these new vari- 
ables, the non-vanishing correlations from Eq. © are 



(OAO) = -cr\ 



(xx) = 



15 ' 



(T- 



(10) 



<AOAO> = o\ , (yy) = . . 
The determinant of the covariance matrix then becomes 



det(M) = 



with y := (t\o\ - cr\ , 



6834375 

and the quadratic form, Eq. (0, turns into 

Q = 31/-^ 15I 2 | 5f | 15(& + &+&) 
2<x 2 2<x 2 2<x 2 2<x 2 

cr 2 (AO) 2 2cr 2 OAO cr 2 2 



(11) 



2y 



2y 



2y 



(12) 



In order to find the number density of potential minima, 
we have to invert the relations given in Eq. (O, considering 
that only the diagonal elements of the tensor £ are non-zero af- 
ter transforming to principal axes. After replacing (fj , % 2 , la) by 
(AO, x, y) and changing the integration boundaries accordingly, 
we integrate only over x and y because the Laplacian of the po- 
tential will become crucial in the following discussion, when an- 
other constraint on AO will be introduced. We can now rewrite 
Eq. © as 



n(0, AO)dOd(AO) = C(M + JV 2 )dOd(AO) 



(13) 



with the integrals 



AO/2 



#1 



J dx J dy|(i 3 -jcy 2 )(A<D-25;)(A<D-3x + y)| 



-AO/2 -3.V-AO 

X |(AO + 3x + }T)| e~ fi , 

AO/2 AO/2 



(14) 



N 2 = j dx j dy|(x 3 - icy 2 ) (AO - 2y) (AO - 3i + y)| 



3JE-AO 

x |(AO + 3i+y)|e^ , 
the normalisation constant 

c- 25 ^ 

and the quadratic form 



Q = 



15x 2 5f o" 2 (AO) 2 2cr 2 OAO o- 2 2 

— + ^2 + — + — + 

o-j o-j y y y 



(15) 



(16) 



(17) 



Equations (1 1 4-1 15b can be integrated analytically, giving iden- 
tical results. The final expression for n(0, AO) is 



n(0,AO) = 



2407r 3 cr 3 Jl5^ 



exp 



(2cr 2 AO + cr 2 o) O 



2y 



x (Fi + F 2 ) , 



(18) 



where F\ and F 2 are functions depending only on the field's 
Laplacian, but not on the field itself, 



= 2cr 2 (5A0 2 - 16cr 2 )exp 
+ cr 2 (l55A0 2 + 32o-^)exp 



(6cr 2 cr 2 -5cr?)A0 2 



2o*y 
(9o- 2 cr 2 -5o-4)A0 2 



8cr 2 



27 



(19) 



F 2 = 5 VlO^AO (AO 2 - 3cr 2 ) exp 



o- 2 A0 2 ^ 



2y 



erf 



V5AO 
2 V2< 



CT 2 J 



+ erf 



' V5AO) 
V2o- 2 J 



(20) 



We point out that Eqs. ( fT8T[20l i are valid in this form only for 
AO > and O < because the underlying integrations over x 
and y were carried out under these restrictions. Both assumptions 
are appropriate; the first because of Poisson's equation, and the 
second because we are only interested in gravitationally bound 
objects whose potentials must be negative. 

For the further evaluation of Eqs. ( fT8T[20| i, we need the first 
three spectral moments of the potential power spectrum, defined 
in Eq. ©. 

3. Linear and non-linear evolution of gravitational 
fluctuations 

The potential power spectrum P^,(k) is related to the density 
power spectrum Ps(k) through Poisson's equation. The power 
spectrum, however, only describes the linear evolution of fluc- 
tuations for which 6 5= 1. Thus, we also need an ansatz for 
their non-linear evolution having higher amplitude. We shall use 
the spherical-collapse model (SCM) to model non-linear effects. 
Along the way, we shall introduce a proper definition of a filter- 
ing scale R. 



3. 1 . Linear power spectrum 

The gravitational potential is related to the density contrast in 
comoving coordinates by Poisson's equation 



AO = 4nGp b a 2 6 , -£ 2 = 4nGp b a 2 6 , 



(21) 



in real and Fourier space, respectively. By the definition of the 
power spectrum, 



{6(k)6*(k')) =: (2nfP(k)5 D (k - k') , 



(22) 



where £ D denotes Dime's delta distribution, and using p b = 
(3// ( 2 f2 m o)/(8^Ga 3 ), the potential power spectrum is related to 
the density power spectrum by 



9 & n HI 



(23) 



Since Ps{k) oc k for k «: ko and Ps(k) oc k for k » ko, where ko 
is the comoving wave number of the perturbation mode entering 
the horizon at matter-radiation equality, we have P<i,(k) oc £~ 3 for 
k ^ ko and P<$(k) oc k~ n for k » ko- 

Due to the steepness of the power spectrum, we have to in- 
troduce a cut-off wave number k m \ n when evaluating the spectral 
moments, 



k 2 dk 
~2n 2 



■P®(k)W 2 R (k) . 



(24) 
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Thus, k m [ n defines a sharp high-pass filter in £-space. It has to 
be chosen properly to filter out large potential modes and there- 
fore also large-scale potential gradients responsible for peculiar 
velocities of collapsed structures. In this way, this filter ensures 
that the gravitational potential of a structure is defined with re- 
spect to the large-scale potential value in its direct vicinity and 
that the constraint of a vanishing potential gradient is fulfilled for 
structures of all sizes. If they moved, they would not be counted 
when searching for potential minima and the number density de- 
rived in that way would be too small. We will discuss later how 
to find the proper fe m j n . 

The evolution of the density power spectrum between the 
scale factors a\ and a 2 is parametrised by the linear growth fac- 
tor D + (a) and the transfer function T(k, a), 



P s (k, a 2 ) = 



D + (a 2 ) 



D + (ai) 



T 2 (k,a 2 )P s (k,ai) 



(25) 



Since the transfer function T(k, a) only changes for redshifts z <: 
100, we do not need to take it into account for the evolution 
of the power spectrum at z < 100. Together with Eq. d23l ). the 
evolution of the potential power spectrum is thus given by 



P®(k, ai) = 



G+(a 2 ) 



G + (fli) 



P®(k, a\) 



(26) 



which defines the potential growth factor G+{a) :- D + (a)/a. 
We normalise G+ such that G+(a = 1) = 1. Since D+(a) = a 
in an EdS universe, the potential growth factor stays constant, 
thus the potential power spectrum does not evolve with time in 
this case. This is not true for a ACDM and an OCDM model, 
yet the variation with time remains small. For cosmologies with 
Q m + 1, Qa + and negligible radiation density, the potential 
growth factor is accurately approximated by 

G + (a) = = 5 -^a)W(a)-D. K {a) 



1 + -Q m (a) 



1 + — O a (a) 
70 



(27) 



(ICarroll et al.lll992h . Figure Q] shows the evolution of G+ with 
redshift. Obviously, the expression "growth factor" is somehow 
misleading because the potential power spectrum's amplitude is 
in fact decreasing with time. 

3.2. Non-linear evolution 

Aiming at the number density of potential minima including 
non-linear evolution, we have to relate the potential from linear 
theory used so far, O], to the potential including non-linear evo- 
lution, <t>„i. We shall use the spherical-collapse model to estimate 
the ratio Oni/^. 

3.2.1. Gravitational potential in the centre of a homogeneous 
overdense sphere 

A spherical and homogeneous overdensity with density contrast 
6 and radius R has a gravitational potential given by Poisson's 
equation, 



r 2 dr \ dr 



: 4nGp0(R - r) 



(28) 
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Fig. 1. Potential growth factor for three different cosmologies. 



where 8 is Heaviside 's step function and p is the density inside 
the sphere acting as the source of the gravitational field. This 
equation holds in physical coordinates, while Eq. ( f2Tb uses co- 
moving coordinates. Integrating twice with the boundary condi- 
tions (<9<3>/<9r) — > for r — > and <$> — > for r — » 00 yields 



O(r) 



-2nGpR 2 \l 

4nG_R? 
3 r 



r 

3R 2 



for r < R , 
else , 



showing that the potential at the centre, <t>o := <t>(r = 0), is 



-2nGpR 



(29) 



(30) 



Since the potential cD(r) oc 1 - r 2 /(3R 2 ) inside a homogeneous 
sphere, it may appear appropriate to construct the low-pass filter 
from the window function 



(31) 



instead of a top-hat. It Fourier transforms to the filter function 



5 3 sin u - u ( 3 + u 2 ) cos u 
W R (k) = -t '- i with u — kR . 



2u 5 



(32) 



Since the window function is more compact in real space than 
the top-hat, its Fourier transform is slightly broader in £-space 
and thus includes more Fourier modes. 

This window function does not (and should not) reproduce 
the potential outside an isolated, homogeneous sphere, which 
drops oc r . This is no surprise because it drops to the poten- 
tial of the homogeneous background within a finite radius, and 
the presence of a background potential signals the breakdown of 
Newtonian gravity in the cosmological context. 

3.2.2. Filtering radius 

The preceding consideration also provides a proper definition 
for the filtering radius of the window function Wr used for the 
calculation of the spectral moments, Eq. (|6]). Poisson's equation 
for the perturbations in physical coordinates reads 



A r <I> = 4nGp 



(33) 
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with p - pbS. Combining Eqs. ( f30l > and d33l l, we see that the 
central potential is Oo = -jA r 0^ 2 . We can use this relation to 
define a filtering radius in physical coordinates by 



Rr 



-20 

aTo 



(34) 



This expression remains valid in comoving coordinates if we re- 
place the Laplacian in physical coordinates by the Laplacian in 
comoving coordinates, 



Rr 



-20 
a 2 A r O 



-20 

A conl O 



(35) 



3.2.3. Spherical collapse model 

We rescale the scale factor a and the radius of the overdense 
region R by their values at turn-around, a ta and R ta , respectively, 
defining the parameters 



x :- 



a R 

"la /Vfl 



(36) 



In addition, we introduce the dimensionless time r and the over- 
density at turn-around I by 



t :- H u t , p ta = |pb,ti 



(37) 



with the Hubble parameter at turn-around H ti . In the follow- 
ing, we will sketch the most important steps to consider cos- 
mologies like EdS, ACDM , and O CDM. Our consideration is 
based on Bartelmann et al. (2006), but simplified because we 
ign ore dynamical dark energ y for now. This extends the work 
by IWang & Steinhardt (1998) towards cosmologies with non- 
vanishing curvature. 

Spherical collapse is then described by the differential equa- 
tions 



+ Q A . ta x 2 + (1 -fi n 



^A,ta) 



1/2 



y 



fira.taf 

2y 2 



(38) 



(39) 



with the matter and dark-energy densities at turn around, D. m jta 
and QA,ta> respectively. Primes denote derivatives with respect to 
the dimensionless time r. Equations ( l38l and ( 1391 can be solved 
numerically using the boundary conditions y — at x = and 
y 1 = at x = 1, meaning that the overdensity starts with zero 
radius and reaches its maximal extent at turn-around. The re- 
quirement y — 1 at turn-around then uniquely determines f . 

The overdensity inside a halo with respect to the background 
density at any time is A = (x/y) 3 %. Assuming that the collapsing 
halo virialises at the collapse redshift z c , which corresponds to 
the normalised scale factor x c when it would ideally collapse to 
zero radius, the virial overdensity is 



1 +Zc 



(40) 



The redshifts at collapse and turn-around, z c and z ta , are related 
by f(zta) = jt(zc). Later, we shall need the overdensity with re- 
spect to the critical density rather than the background density, 
A v = O m (z c )A v . The ratio between the radii at turn-around and 
at collapse is 2 in an EdS universe because of virialisation, but 
slightly different in more general cosmologies. Due to an ad- 
ditional contribution of dark energy to the potential, denoted by 



(E^), the virial theorem is modified to (Ziyn) = -j{Epot) + {Ede)- 
This l eads to the following app roximation for the ratio of the two 
radii dWang & Steinhardll998l) : 



R v 1 2 r l v 

«ta 2 + 77t - \t] v 



T]y = 



2 Qa(zc) (l+z 



£ Q m (z c ) \ 1 + z t 



2 OAfaa) 
I ftm(Zta) 



(41) 



As the density perturbations evolve linearly at very early times, 
the density contrast expected from linear theory at collapse, 
6 c (x c ), is simply given by 



«5 c (x ) = hm(^^[A v (x)-l] 



(42) 



In an EdS universe, 6 C = 1.686 and A v = 178, both independent 
of the collapse redshift. 

3.3. Ratio of linearly and non-linearly evolved potential 

We shall use the potential in the centre of a spherical and homo- 
geneous overdensity, Eq. ( l30t . to relate the linear and non-linear 
potential depths. During the following calculation, we will de- 
note quantities at an initial scale factor a; with the subscript T 
and quantities at the collapse scale factor a c with the subscript 
'c'. 

The potentials at the initial time and at collapse are 



Oi = -InGpfiiyl . ~- 
respectively. Their ratio is 

£c _ <5v /yc\ 3 /xA 3 yj_ 
<J>i ~ Si \x c ) [yij y c ' 



-2nGp c Rlyl 



(43) 



(44) 



where 6 V = A v - 1. We have used here that the densities at 
the initial time and at collapse time are p; = pi,o^~ 3 ^i and 
Pc = Pbo<2c 3 £> v , respectively. Using Eq. d40] >. we can write 
(yc/xc) 3 = ^/A v and y- t « £ 1/3 Xi since A; * 1 for early times. 
This yields 



s,t 



1/3 



1/3 



A v y c Si y c Si 



(45) 



We have neglected the difference between S v and A v in the last 
step, which is a good approximation since A v = <9(10 2 ). The 
ratio Si/xi =: C is given by 



? l/3 h 



A.ta 



1 — £2 m , ta — £Ia 

fim.ta 



(46) 



<f^m,ta / 

(Bart elmann et alj|2006i) . Since Eqs. ( 145146b describe the non- 
linear evolution between scale factors aj and a c , the linear and 
non-linear potential depths are related by 



0„i _ O c Oj 
Oi O; Oi c 



<r 1/3 1 G+joi) 

y c CG+(a c ) ' 



(47) 



where Oi ;C is the potential evolved linearly from «i to a c . For an 
EdS universe, we have O n i/Oi = -j, independent of the collapse 
redshift because the potential growth factor G+(a) is constant 
and y c = j, but for more general cosmological models, it de- 
pends on the collapse redshift. 
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Fig. 2. The ratio of the potential depth that is expected at collapse 
redshift z c from linear theory, Oi, and the one inferred from the 
SCM, 0>„i. 



We show O n i/Oi for three different cosmological models in 
Fig. [2] We arbitrarily choose the initial scale factor to be five 
times the scale factor at matter-radiation equality because this 
is early enough for the potential not to have begun developing 
since we are in the matter dominated era, and late enough for 
ignoring any evolution of the transfer function T(k, a) with time. 
While the ratio Eq. (|47| > reaches the constant expected in an 
EdS model relatively quickly in the ACDM case as the redshift 
increases, its evolution is much slower for OCDM. We point out 
that the difference between the linear and the non-linear potential 
evolution in the centre is by far not as large as for the density 
contrast, where it is of order 10 3 4 . 



3.4. Counting only collapsed structures 

Since we are only interested in collapsed structures forming ha- 
los, we have to include this additional pro perty in our calcula- 
tions. Following Press & Schechter (1974), this can be done by 
only taking structures into account whose linear density contrast 
exceeds the critical value <5 C . In our case, this translates to the 
criterion that the Laplacian AO exceed a certain threshold A<J> C , 
which is, according to Poisson's equation, given by 



«5c(a) 



AO c (a) = -H£Q m o 



(48) 



In that way, we only count structures whose gravitational poten- 
tial is "curved enough", so that the number density of potential 
minima belonging to collapsed structures is finally 



n(O) 



■/ 

AU>, 



d(AO)n(0, A<D) 



(49) 



with «(€>, AO) given by Eqs. 



4. Construction of the X-ray temperature function 

Combining the results of the previous sections, we are now able 
to derive a function describing the differential number density of 
structures as a function of their X-ray temperature based on the 
statistics of minima in the cosmic gravitational potential field. In 
order to count them properly, a high-pass filter is introduced that 



removes disturbing large modes. Assuming virial equilibrium, it 
is possible to relate the potential depths to X-ray temperatures 
using the virial theorem. 

4. 1 . Virial theorem 

We can relate the potential depth O n i to the X-ray tempera- 
ture of a cluster via the virial theorem for Newtonian gravity, 
which relates the ensemble-averaged kinetic and potential ener- 
gies, (E kin ) and <£ pot >, by (E^) = -±<£p 0t >. The kinetic energy 
is connected to a temperature T by (£kin) = jk^T, where kn is 
Boltzmann's constant. The potential energy is mO with a proper 
mass m. Assuming that the intracluster medium is fully ionised, 
we know that m = ymi p with fi - 0.59. Particles near the cluster 
core feel the potential O(r) w Oo, thus the virial theorem read^] 



(50) 



Since Oo can be regarded as the non-linear depth of a po- 
tential minimum, we can replace O n i by <t>o = -3IcbT /(pm p ). 
Given a particular X-ray temperature, we calculate the corre- 
sponding linear potential depth by relating the temperature to 
the non-linear potential using Eq. 150) . and the latter to the lin- 
ear potential using Eq. (147) . arriving at 



O, 



3y c C G+(z) 
"fi/3 jump G + ( Zi )' 



(51) 



4.2. Evaluation steps 



The X-ray temperature function is implicitly determined by 
Eqs. ( fT8l420b . Assuming a temperature T, the linear potential 
depth <t>i is found from Eq. ( Bit . Then, Eq. ( TT~8T > needs to be inte- 
grated over AO from AO c to infinity. Since the smoothing radius 
(l35T l depends on the integration variable AO, each step in the 
integration requires updating the spectral moments, cf. Eq. (0. 

In practice, the temperature interval for which the X-ray tem- 
perature function is to be calculated is divided into a reasonable 
number of bins large enough that the shape of the temperature 
function can be inferred from interpolating between them. At 
first, each X-ray temperature T, corresponding to one particu- 
lar bin has to be related to a linear potential depth Oi , using 
Eq. ( Bit . We also have to take care of the Jacobian determi- 
nant of the transformation because n(T)dT must be replaced by 
«(Oi)(dr/dOi)dOi. Having related both quantities, it is possible 
to evaluate Eqs. (fT8ti20b . 

In order to calculate the spectral moments <tq, ct\, and o"2, 
we also have to choose a reasonable amount of bins for the 
Laplacian of the potential, which we will denote with AOy, since 
the spectral moments are functions of both Oi and AO through 
the filtering radius entering via the Fourier transform of the win- 
dow function, Eq. ©. Additionally, the cut-off wave vector k m [ n 
defining a sharp high-pass filter in A:-space is also a function of 
both quantities. A detailed discussion how to find the proper k m j n 
for a given temperature and Laplacian of the gravitational poten- 
tial is presented in Sect. 14.31 Thus, we have to evaluate both the 



1 Strictly speaking, this is only valid for a universe with C1 A = 0. 
Due to the presence of dark energy, an additional potential arises 
(Wang& Steinhardt 1998) whose contribution is small (e.g. the virial 
radius only changes about few percents, see Sect. l3.2.3t . so that we ne- 
glect it in our further calculations. Whether dark energy contributes at 
all to the pro cess of virialisation is st ill an open question, see e.g. the 
discussion by Bartelmann et al. (2006). 
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filtering radius and k m [ n for each pair (Oi AO^) using Eq. < f33T > 
and starting for a given Oy with AOo = AO c . 

Inserting the spectral moments into Eqs. ( [T8l - l20l i yields 
the number density of minima per potential interval dOi and 
per interval of the Laplacian d(AO) for the specific parame- 
ter pair (O u , AO c ). Since A<D C = 0(1Q 4 km 2 s~ 2 Mpc~ 2 h 2 ), 
we choose the next step to be A<E>i = AO c + c)(AO) with 
(5(AO) = 10 2 km 2 s~ 2 Mpc~ 2 h 2 and calculate now R and the 
spectral moments for the pair (Oi.,, AOi). We continue with 
AO, = A4> 7 -_i + £(AO) for a given Oy until we fulfill the fol- 
lowing convergence criterion. 

We approximate the integral in Eq. d49l ) by the trapezium 
rule. Hence, for a specific potential value Oi the number density 
n(<E>i,,) is evaluated numerically as 

«(a>i,i) * 2 [" (01 -" Aa> >- i} + " ( °''" A ®M ■ (52) 

This summation is stopped at an index N chosen as the first index 
for which the relative contribution to n(Oi,,) is smaller than 1CT 6 . 
This is a proper break condition since AO) tends rapidly 
towards zero for AO — > oo. 

Having evaluated n(Ti) for each i, we have an appropriate 
approximation of the X-ray temperature function's shape in the 
chosen temperature interval. 



4.3. High-pass filtering 

As mentioned before, we need a proper high-pass filter in order 
to remove disturbing large-scale potential modes since we want 
the potential associated with a collapsed structure to be defined 
with respect to the large-scale potential level in its direct vicinity, 
and we want the structure to have no peculiar motion so that the 
constraint rj — is applicable. The natural filter choice is a sharp 
cut-off in £-space since this will effectively remove both large- 
scale potential modes and potential gradients. 

Let again R be the top-hat filter radius of Eq. (|6]l and R^ p the 
filter radius related to a sharp cut-off wave number &h p in fc-space 
by &hp = 27r/7?h p . The ratio of both quantities defines a, 



fl hp = aR 



(53) 



Since R and R^ p define a low- and a high-pass filter, respectively, 
one can expect a > 1 . 

It is now a legitimate question how a should be chosen for 
calculating the correct number density of objects having a par- 
ticular potential depth. We argue that a natural choice exists. 
In Fig. [3] we plot h(T, AO) as a function of a. For each pair 
(Tj, AOjt), the number density peaks at some value a max , which 
is a function of both the temperature and the Laplacian of the po- 
tential. It increases with both increasing Laplacian and increas- 
ing temperature. For low temperatures, the maximum is rather 
sharp, but broadens as the temperature increases. This is most 
pronounced for the EdS universe. This behaviour can be under- 
stood as follows. 

On the one hand, decreasing the radius of the high-pass filter 
starting from a value much larger than the typical size of the ob- 
ject excludes more and more modes on decreasing scales. In this 
way, potential gradients are removed which would cause a non- 
zero peculiar velocity and thus a deviation from the constraint 
r\ — 0. Since more and more objects are put to rest, the number 
density of objects with vanishing potential gradient increases. 
Figure [3] shows that the increase of the number density is less 
steep for large (or hot) than for small (or cool) objects. This re- 
flects the fact that that massive, hot objects, e.g. with T — 10 keV, 



2A<1> 




10 12 14 16 18 20 



2A<1> C 
3A<1> 

4A[>: 




10 12 14 16 18 20 




14 16 18 20 




Fig. 3. Differential number density h(T, AO) as function of the 
parameter a for three different cosmologies and two tempera- 
tures at z = 0. Upper row: ACDM. Central row: OCDM. Lower 
row: EdS. Left panels: T — 1 keV. Right panels: T = 10 keV. 
The values for AO were chosen to be multiples of the critical 
Laplacian AO c and are given in the plots. For T — 10 keV we 
only plot the number density for AO = AO c since it is a very 
steep function of AO and it is too small for larger multiples of 
AO r to be seen. 



are more likely located at potential minima, thus removing large- 
scale modes has less effect on their number counts. 

On the other hand, increasing the radius of the high-pass fil- 
ter starting from a value much smaller than the typical size of the 
relevant object adds more and more modes. While the window 
between low- and high-pass filtering is too small, modes rele- 
vant for the structures considered are filtered out and the halo 
number density remains approximately zero. Once modes are 
included that compose the structures, the halo number density 
steeply rises as a is increased. 

At a certain a max , the number density of objects reaches a 
maximum where both effects are balanced. Then, all modes rel- 
evant for structures of size R = -y/-20,/AO J - are included, but 
larger modes are excluded which would create a non-vanishing 
potential gradient. Thus, a max can be used to define k m i n for each 
pair (T,-, AOj) individually by 



2n 

CUmaxR 



(54) 



This definition of k m [ n has to be used when evaluating the spec- 
tral moments with Eq. ( f24b . 

In appendix [A] we present an alternative way to define a 
physically reasonable cut-off wave vector for the evaluation of 
the spectral moments. It turns out, however, that the number den- 
sity of objects with a low X-ray temperature function is highly 
underestimated in this way. 
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4.4. Inferring the temperature function from the 
Press-Schechter mass function 

In order to compare the X-ray temperature function that we have 
derived from the statistics of gravitational potential perturbations 
to the classical Press-Schechter theory, we need a proper, albeit 
idealised, and consistent mass-temperature relation. Since we 
used the virial theorem, Eq. (l50l >. to relate the temperature to the 
potential, we will start at the same point to relate the temperature 
to a mass. Note that this has nothing to do with an assumption 
on real clusters, but merely serves the purpose of a theoretical 
cross-comparison between the mass-based Press-Schechter ap- 
proach and our direct derivation of the temperature function. 

We saw earlier that for a spherical and homogeneous over- 
density the potential depth in the centre is <t>o = -2nGpR 2 , 
where p = 6 v pb is the constant density inside the perturba- 
tion. Nonetheless, we shall replace 5 V by the virial overdensity 
A v , which is a good approximation because 6 V = A v - 1 and 
A v = O(l0 2 ). 

The mass of the overdensity is M — ^npR 3 , where p is the 
total density inside the sphere. It is related to the background 
density by p — A v pb- According to the previous statements, we 
can identify p and p, thus p as p. Combining the equations for 
the potential and the mass and using the virial theorem yields the 
temperature-mass relation 

k B T = Gpm p M 2 ' 3 . (55) 

The density inside the cluster p is related to the background 
density by p(z) = p cr o^moA v (z)(l + z) 3 . In our further calcula- 
tions, however, we will use the virial overdensity A v , which was 
defined with respect to the critical density instead of the back- 
ground density. This yields 

P(z) = PcM -^-K(z)(l + z) 3 . (56) 

Inserting the previous equation into Eq. d55l l. normalising to 
M = 10 15 M Q and A v = 178, and combining all other quanti- 
ties including Hq = lOO/ikms -1 Mpc~' into one normalisation 
factor, we find 

/ Q m0 A v (z)\ 1/3 / M \ 2/3 

This is almost the same relation as given bv lEkeetal.l(ll996h . 
The only difference is the normalisation factor. While we used 
a spherical and homogeneous overdensity for the calculation in 
order to be consistent with the derivation of the potential func- 
tion, they derived a mass temperature-relation for an isothermal 
sphere, which results in a different normalisation. 

We shall use this relation to convert the classical Press- 
Schechter mass function into a temperature function. Again, we 
must not forget to account for the Jacobian when transform- 
ing from mass to temperature. In contrast to the transformation 
of the potential to the temperature, the Jacobian of the mass- 
temperature transformation depends on temperature since both 
quantities are non-linearly related. 

Finally, we also use Eq. ( [57]) to conve rt the mass functions 
derived by ISheth & Tormenl d 19991) and iJenkins et all d200ll) 
to X-ray temperature functions knowing that these are fits to 
numerical simulations including effects of ellipsoidal collapse, 
while the derivation of our mass-temperature relation is based on 
spherical collapse. Yet, it should yield qualitative information on 
the importance of ellipsoidal collapse in our potential approach, 
since structures in the potential field are smoother than structures 
in the density field. 



5. Results 

In this section, we present the X-ray temperature functions for 
the ACDM, OCDM, and EdS cosmologies calculated from the 
statistics of gravitational potential perturbations and compare 
them to temperature functions derived from three well-known 
mass functions. The results are shown in Fig. |4j where we com- 
pare the X-ray temperature function inferred from the potential 
using Eq. ( BTT l to the X-ray temperature function calculated from 
the Press-Schechter, the Sheth-Tormen and the Jenkins et al. 
mass functions using Eq. d57l ) for our three cosmological models 
and three redshifts. 

The X-ray temperature function inferred from the statistics 
of fluctuations in the gravitational potential matches the temper- 
ature function inferred from the Press-Schechter mass function 
quite well for all three cosmological models and redshifts shown. 
As expected, the filter modelled after the internal potential pro- 
file of a homogeneous sphere (filter 2) yields slightly larger num- 
ber densities especially at the high-temperature end because it is 
wider in fc-space than the top-hat filter (filter 1). The temperature 
function based on our novel approach is in good agreement with 
the classical Press-Schechter approach with an only slightly dif- 
ferent amplitude depending on redshift. These differences, how- 
ever, may be irrelevant because of the idealising assumptions 
entering both approaches. 

A comparison with temperature functions inferred from 
mass functions including elliptical collapse like the Sheth- 
Tormen mass function shows that the deviations increase sub- 
stantially, especially for high redshifts. This suggests that ellip- 
soidal collapse should also explicitly be included in the poten- 
tial approach. However, we emphasise again that this compari- 
son is only qualitative since the mass-temperature relation used 
is based on spherical and homogeneous objects so that the ac- 
tual difference between the different approaches may be even 
smaller. 

6. Discussion and Conclusions 

In this paper, we developed a novel approach to a theoretical 
derivation of an X-ray temperature function that does not de- 
pend on global, unobservable cluster quantities, but merely on 
the depth of the gravitational potential, which is a locally de- 
fined quantity directly related to observables. Using the sta- 
tistical properties of a Gaussian random field, we were able 
to derive a distribution for the minima in the cosmic gravita- 
tional potential. Counting only those potential minima which 
are "curved enough", in the sense that their Laplacian exceeds 
a critical value, it was possible to semi-analytically compute a 
potential function belonging to collapsed structures. The crit- 
ical Laplacian needed to distinguish collapsed structures from 
non-collapsed structures could be related to the critical density 
contrast 6 C of the spherical-collapse model, which also plays an 
important role in the Press-Schechter formalism. 

We also managed to calculate the influence of non-linear 
structure formation on the potential by referring again to the 
spherical-collapse model, which allows the ratio of the linearly 
and non-linearly evolved potential depths to be computed. Both 
the linear and the non-linear evolution of the potential are much 
slower compared to the evolution of the matter density. The ratio 
between the non-linearly and linearly evolved potential depths, 
for example, is ~ 30 times smaller than the ratio A V /6 C . 

One crucial ingredient to our approach is a proper high-pass 
filter that removes large-scale potential modes and gradients re- 
sulting from them. We find that choosing the filter scale such as 
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Fig. 4. Comparison of the X-ray temperature functions derived from the statistics of gravitational potential perturbations (<J> n i) and 
the mass functions by Press & Schechter (PS), Sheth & Tormen (ST), and Jenkins et al. (J) for three different cosmologies and three 
redshifts. The filters 1 and 2 are defined by Eqs. d6} and d32b . respectively. Upper row: ACDM. Central row: OCDM. Lower row: 
EdS. Left panels: z = 0. Central panels: z = 1. Right panels: z = 2. We want to emphasise that the y-axes are not scaled equally in 
order to compare the temperature functions more appropriately. 



to maximise the number density of objects yields good agree- 
ment between the Press-Schechter approach based on an ide- 
alised mass-temperature relation, and our direct derivation of the 
temperature function. This criterion thus provides the foundation 
for using number counts of galaxy clusters in cosmology with- 
out invoking global quantities like cluster masses. Although we 
believe that our method is essentially ready to be tested with sim- 
ulations before applying it to observational data, it leaves room 
for improvements in the theoretical description, as follows. 

Based on this work, it should be examined if elliptical col- 
lapse can improve the results compared to fully non-linear N- 
body simulations than spherical collapse. First, we need to find 
out if the forma l ism f or th e elliptical collaps e developed by 
ISheth & Tormenl d 19991) and lSheth etafl d200ll) can be adapted 
for the potential calculation or if another approach must be 
found. Second, it will be interesting to see whether this will 
influence the results significantly, as suggested by Fig. |U This 
is not as obvious as for density perturbations because potential 
fluctuations are much less asymmetric, and their non-linear evo- 
lution is much less pronounced than for density perturbations. 

We calculated the X-ray temperature function from the po- 
tential function simply by applying the virial theorem relating 



the kinetic energy to the potential energy. It is straightforward to 
relate a temperature to a potential depth in this way. But there are 
two issues that should be examined in more detail. The first point 
is the influence of dark energy on the virialisation process. We 
have already mentioned that it might not affect virialisation at all, 
but if it does, it should be relatively small since the variation of 
the virial radius due to the additional potential originating from 
the dark energy is at the level of few per cent. A better under- 
standing could strengthen the assumption that it can be neglected 
when relating the potential to a temperature for models includ- 
ing dark energy like ACDM. The second point is that we used 
the potential at the minimum of the potential well in the virial 
theorem. But since the virial theorem has to be applied to an 
averaged potential (<D) instead of the potential in the minimum 
<t>o, this should result in an additional correction factor when us- 
ing <t>o instead. A closer examination of this factor, especially 
its magnitude and its dependence on the potential depth, needs 
to be carried out. Besides, real clusters have lower temperatures 
in their centres due to cooling, implying that the application to 
observations requires calibration. 

On the whole, our results are very promising and suggest 
to continue following this approach, which should allow a di- 
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rect comparison of the cluster population with cosmological pre- 
dictions based on observable, local cluster quantities. Although 
we only presented results for ACDM, OCDM, and EdS, it is 
straightforward to extend our computation of the X-ray temper- 
ature function to more elaborate cosmological models including 
e.g. quintessence or early dark energy. 
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Appendix A: An alternative way to determine the 
cut-off wave number 

In this appendix, we present an alternative approach to determine a physically 
reasonable definition of k m \ n . Although it does not give the correct number den- 
sity for smaller objects with a low X-ray temperature, the definition presented 
here may become important for future work based on the potential perturbation 
approach. 
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Fig. A.l. Comparison of the X-ray temperature function derived 
from the statistics of gravitational potential perturbations (<t> n i) 
using the alternative definition of k m m from Eq. dA.lt with the 
classical Press-Schechter approach (PS) for three cosmologies. 
Upper row: ACDM. Central row: OCDM. Lower row: EdS. Left 
panels: z — 0. Right panels: z — 2. 



A.1. Definition 

An alternative appropriate choice for k m i n could be the redshift-dependent parti- 
cle horizon r^ OI , taking into account that light could have travelled only a finite 
comoving distance between the Big Bang and redshift z. Consequently, we must 
only consider modes of the gravitational potential already inside the horizon. 
Thus, 



nior(z) C 



lim f — 

Ji->oJ a '- 



da' 



E(fl') 



(A.l) 



where E(a') is the expansion rate of the Universe evaluated at the scale factor a' . 

In principle, signal retardation should also be taken into account. 
Considering an arbitrary point xq at time to, only modes lying inside its past 
light cone can have influenced it. Retardation has to be included "by hand" be- 
cause we are using Newtonian gravity. This gives rise to an additional factor 
when calculating the power spectrum's amplitude because we must evaluate the 
amplitude of a mode with wave length A not at time to, but at the earlier time 
t - At = t - A/(2c). 

We can compute the corresponding scale factors at to - At as follows. We 
must have A = 2n/k = 2D com (z, zk), where z and zt are the redshifts correspond- 
ing to times to and to - At, respectively, and D com (z, Zk) is the comoving distance 
between both redshifts. Thus, we have to find a scale factor ak = 1/(1 + z,k) for 
each mode k such that 



nH 



D com (a,a k ) 



/ 

Vi k {- k ) 



da' 



a' 2 E(a') 



(A.2) 



holds. This is consistent with Eq. jA.lt because the wave number k approaches 
^min f° r a k —* 0. The influence on the power spectrum's amplitude results in an 
additional factor G\(cik)IG\(ci) entering Eq. )26t . Since P$ does not evolve with 
time in an EdS universe, it has no effect in this case. Additionally, it turns out that 
its contribution is quite small for both the ACDM and the OCDM model, where 
it only affects the power spectrum's amplitude for a relatively small amount of 
wave numbers. Hence, it only negligibly affects the computation of the spectral 
moments and can usually be ignored. 



A.2. Results for the X-ray temperature function 

In Fig. IA.ll we present the results for the X-ray temperature function from the 
statistics of gravitational potential perturbations using the alternative definition 
of 

^min an d compare it to the Press-Schechter approach for three different cos- 
mologies and two redshifts for temperatures between 1 and 30 keV. We can see 
clearly that both functions match quite well for very high temperatures, espe- 
cially in the EdS case. For all three models, the temperature function derived 
from the gravitational potential is much flatter for low temperatures than the 
temperature function inferred from the Press-Schechter mass function so that the 
number density of objects with temperatures of about 1 keV is too low by a factor 
of more than 100. 

This discrepancy can be explained considering that we expect many more 
density maxima than potential minima for the same volume of space due to the 
following reason. The Gaussian random field of potential perturbations is much 
smoother and has much more power on large scales than the corresponding field 
for the density contrast due to the the steepness of the potential power spec- 
trum. Only large structures that have a high density contrast also correspond to 
a potential minimum, smaller structures only correspond to a maximum in the 
potential's Laplacian but not to a minimum in the potential itself. Since they are 
not located at a minimum of the potential, they have a non-vanishing potential 
gradient which corresponds to a non-zero peculiar velocity. 

The definition presented in Eq. (54) does not involve these problems due 
to the fact that large-scale potential gradients are removed and therefore, the 
condition i] = is also applicable for structures with a low X-ray temperature. 



